Fast projectile stopping power of quantal multi-component strongly coupled plasmas 
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The Bethe-Larkin formula for the fast projectile stopping power is extended to multi-component 
plasmas. The results are to contribute to the correct interpretation of the experimental data, 
which could permit to test the existing and future models of thermodynamic, static, and dynamic 
characteristics of strongly coupled Coulomb systems. 
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Stopping power is a characteristic of primary interest 
for different areas of physics such as nuclear physics, con- 
densed matter physics and plasma physics, as it arises 
when studying the interaction of charged particles with 
matter. In 1930 Bethe derived his seminal formula for 
the fast projectile energy losses assuming that the atoms 
of the medium behave as quantum-mechanical oscillators 
[H|. Later, Larkin Q showed that when fast ions perme- 
ate an electron gas, an analogous formula is applicable, 
but with the mean excitation frequency replaced by the 
plasma frequency lu p : 
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where In A = \n2mv 2 / fuv p is the quantal Coulomb log- 
arithm, Z p e and v stand for the charge and velocity of 
the projectile, p is the target density, A the mass of the 
target atoms, N the Avogadro number, vf the electron 

1 /2 

Fermi velocity, and uj v — (Aime 2 /m) , m and n being 
the electron mass and density. This formula is usually 
employed to determine experimentally n in a charged 
particle system. Particularly, its applicability seems to 
be more promising in the field of plasma physics [H, 0, HJ 
for two reasons: first, in an ionized medium the energy 
loss is mainly caused by the free electrons, leading to 
an enhancement of the stopping power compared to the 
cold target [1, 0, @ ; secondly, this technique appears as 
the only suitable candidate for the diagnosis of hot and 
dense (n > 10 19 cto~ 3 ) plasmas, because most of the 
other methods fail under these conditions [5j. 

Usually, it is believed that the electronic subsystem 
of a plasma provides the main contribution to the stop- 
ping power process, especially for fast projectiles. Our 
first aim in this Letter is to show that in a multi- 
component completely ionized hydrogen plasma with a 
weakly damped Langmuir mode of dispersion uil (fc) , the 
plasma frequency in the Coulomb logarithm of JT]) should 
be substituted by the long-wavelength limiting value of 
Lu L (k), uj l (0) = ujpy/l + H with H = h ei (0)/3 = 
(<?ez (0) — 1) /3, g e i (r) being the electron- ion radial distri- 
bution function. The generalization to partially ionized 



plasmas or plasmas with complex ions and more species is 
straightforward. This correction may have further prac- 
tical implications, in particular, after the experiments 
reported in Ref. [1] where it was possible to measure 
separately the enhancement of the stopping power of fast 
ions due to the increase in the Coulomb logarithm, In A. 
Thus, this method will permit to probe directly strong 
coupling effects which are relevant to plasmas within the 
high density energy regime. This includes plasmas aris- 
ing in astrophysics and space science, planetary interiors, 
inertial confinement fusion, matter under extreme condi- 
tions, metals and condensed matter plasmas. 

Leaving the ionization losses aside, for calculating the 
stopping power for a fast projectile passing through a 
Coulomb fluid we will adopt the polarizational picture, 
which becomes more accurate as the kinetic energy of the 
projectile increases. In 1954 Lindhard obtained an ex- 
pression relating the polarizational stopping power with 
the medium (longitudinal) dielectric function Q. This 
expression can be generalized further by applying the 
Fermi golden rule to obtain @, H, : 
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a± (fc) = ±kv + hk 2 /2M, where M is the mass of the 
projectile (here we will work with heavy-ion projectiles, 
M ^> m), and tib = (1 — exp (— phw))" 1 , fi~ x being the 
temperature in energy units. In addition, unmagnetized 
Coulomb fluids are considered and, hence, the dielectric 
function effectively depends only on the wavevector mod- 
ulus. Expression J2} is valid only if the interaction be- 
tween the projectile and the plasma is so weak that it 
can be treated as a linear effect and no relativistic ef- 
fects need to be taken into account, i.e., when the energy 
lost by a projectile is much less than its kinetic energy, 
which, in turn, is assumed to be much smaller than its 
rest energy |26j |. 

The literature on the polarizational stopping power is 
very extensive. The problem has been analyzed within 
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the random-phase approximation (RPA) [t| and beyond, 
introducing an analytic formula for the local field cor- 
rection (LFC) factor [lfj. In addition there are also 



nonlinear polarization effects which are beyond the 
scope of this work. Whereas we assume that the cou- 
pling between the projectile and the target plasma can 
be treated perturbatively, we do not impose any restric- 
tion on the value of the coupling parameter, T = (3e 2 /a 
(a = (47m/3) -1 ' 3 being the Wigner-Seitz radius), with 
the proviso that the latter remains in the liquid phase 
[2^ . As said before, here we will focus on a completely 
ionized strongly coupled hydrogen plasma. The model- 
ing of the dielectric properties of this kind of plasmas 
constitutes a difficult problem, because its characteris- 
tic lengths, i.e., Wigner-Seitz radius and Debye radius, 
A_d = (47rne 2 /3) -1 / 2 , are of the same order of magnitude 
(in a strongly coupled plasma T = a 2 /3X 2 D > 1, what 
makes mean field theories, such as the RPA, and pertur- 
bative treatments no longer valid) and, at the same time, 
its electronic subsystem is degenerate. 

The framework.- Our dielectric formalism is based 
on the method of moments (13 , [Uj], which allows to 
determine the dielectric function t(k,uj) from the first 
known frequency moments or sum rules. The sum rules 
we employ are actually the power frequency moments 



of the loss function (LF) C (k, ui) 
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defined as C v (k) = ir^ 1 f^uj" £ (k,uj) duj, v — 0,1,... 

Due to the parity of the LF, all odd-order fre- 
quency moments vanish. The even-order frequency 
moments are determined by the static characteristics 
of the system. After a straightforward calculation 
one obtains [H [H, Gi|: C (k) = (1 - e^^O)), 



K(k) 
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H), 



C 2 {k) = uj 2 p , and C A {k) = 

with K(k) = ((v 2 ) k 2 + h 2 k 4 / (2m) 2 J /uj 2 , (vf) being 

the average squared characteristic velocity of the plasma 
electrons. The last two terms in C4 can be expressed in 
terms of the partial structure factors S a b (k) , a, b = e,i: 
(2<K 2 n)- 1 f™p 2 (S ee (j>)-l)f(p,k)dp, 
) f^°p 2 S e i (p) dp, where we have 
5/12 - p 2 /(4k 2 ) + 



U(k) 

H = (67 
introduced / (p, k) 



(k 2 -p 2 Y\n\(p + k)/(p-k)\/(8pk 3 ). 

The Nevanlinna formula of the theory of moments ex- 
presses the dielectric function which satisfies the known 
sum rules {C 2l/ } 2 =0 (l2l. Ila. Il6j|: 
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where uj\ = uj\ (k) = C 2 /Co, uj 2 — UJ% (&) — C4/C2, in 
terms of a function q = q(k,z), which is analytic in the 
upper complex half-plane Im z > and has there a pos- 
itive imaginary part. It must also satisfy the limiting 
condition: (q(k, z)j z) — > as z — * 00 for Im z > 0. 
In an electron liquid this Nevanlinna parameter func- 
tion plays the role of the dynamic LFC G(k,uj). In 



particular, the Ichimaru visco-elastic model expression 
for G(k,uj) is equivalent to the Nevanlinna function ap- 
proximated as i/r m , r m being the effective relaxation 
time of the Ichimaru model [171 ] . In a multi-component 
system the Nevanlinna parameter function stands for 
the species' dynamic LFC's. In general, we do not 
have enough phenomenological conditions to determine 
that function q{k,uj) which would lead to the exact ex- 
pression for the LF. One might benefit from the Perel' 
- Eliashberg (PE) (3| high-frequency asymptotic form 
0, Ime^X^)- 1 ) ~ (4/3) 1/4 r s 3/4 /3K/c;) 9/2 , 
where r s = ame 2 /H 2 is the Brueckner parameter. 

The corrected Bethe-Larkin formula.- Let us choose 
a model function q satisfying the conditions mentioned 
after the Nevanlinna formula ([3|) that would permit to 
treat the stopping power calculation analytically. If we 
put simply q(k,u>) — i0 + , then we get the following par- 
ticular solution of the moment problem: 



J ±8 (uj) + [6 (uj - uj 2 ) + 5 {uj + oj 2 )\ , 
Zuj 2 



C{k,uj) 
ttC (k) 

(4) 

Physically, Eq. (j4|) describes an undamped collective ex- 
citation mode (Feynman approximation) at uj 2 with an 
additional central peak accounting for hydrodynamic dif- 
fusional processes The applicability of this expres- 
sion is justified provided that the damping of the collec- 
tive excitation is small enough, making this mode to act 
as the main energy transfer channel. Thus we can dis- 
regard the details of the rest of the excitation spectrum. 
If we introduce expression ([4]) into the Lindhard formula 
([2|), it immediately reduces to: 
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where the "cut-off" wavenumbers k\ and k 2 are such 
that the inequality < uj 2 {k) < kv is satisfied with 
v/vp — > 00 and u> 2 (k) understood as the plasma Lang- 
muir mode dispersion law ujL{k). For a weakly cou- 
pled plasma the RPA dispersion law is valid which ne- 
glects the correlational contributions to uj^k): Wi(fc) = 



(lu 2 + (v 2 e ) k 2 + h 2 k 4 / (2m)' 
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Then, if v is asymp- 
totically large, we have ki — uj p /v, k 2 = 2mv/h, and we 
recover the Bethe-Larkin (BL) result [l], Q • Notice that 
in the above-mentioned inequality for uj 2 , we have pre- 
sumed that kv ^> Hk 2 /2M, which is equivalent to disre- 
gard, at most, terms of the order of m/M. Similar terms 
were omitted in the above expressions for the moments 
C 2 and C4, as well. 

To take into account all Coulomb and exchange in- 
teractions in the system analytically, we might use for 
the electron-electron contribution U (k) its long- and 
short-range asymptotic forms, U (k — > 0) ~ —v 2 e k 2 /oj 2 , 
U(k^ 00) ~ -/i ee (0)/3, where v 2 ee = -4E ee /(15nm) 
is defined by the plasma electron-electron interaction 
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energy density E ee of the plasma [§] , h ee (0) be- 
ing equal to the previous expression for U (k), but 
with the function f(p, k) replaced by unity. If 
we interpolate the plasma mode dispersion law as 

uj L (k) = + H) +wk 2 + h 2 k i / (2m) 2 ) V2 with 

w = 2{v 2 ) ~ v ee> then the "cut-off" wavenumber k± 
is modified as k[ — uj' p /v, with uj' p = uj p yl + H, for 
v/vp oo, so that the fast projectile stopping power 
becomes: 



dE 

dx «»«f 
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Here, the correction H stems from the electron-ion corre- 
lation contribution to the moment C±(k) and is also the 
one responsible for the upshift in the value of the Lang- 
muir frequency predicted in the long-wavelength limit 
for an electron-ion plasma with an undamped collective 
mode. Although the accurate calculation of H under 
realistic conditions is a difficult task 2(1, 21], it is possi- 
ble to find a simplified analytic expression based on the 
temperature Green's function technique by a regularized 
summation over the Matsubara frequencies [12 1, yielding 
H = (4/3)r s VT/(2 A /rJ + IV6) (see also RefTj^ for an 
alternative approach based on a nontrivial renormaliza- 
tion via pair-correlations in liquid metals). Whereas in 
a weakly coupled plasma, r < 1, this correction is neg- 
ligible, in a strongly coupled Coulomb system it could 
be possible to retrieve directly H (or g e i(0)) by fitting 
Eq. ([6]) to some experimental data. For instance, if we 
take ff e i(0) = 10 [iH and In A = 14 Q, then the stop- 
ping power obtained by the BL formula gets modified by 
~ 5%, which indicates to what extent the experimental 
accuracy needs to be improved. 

The damped collective mode.- The collective mode is 
expected to be damped [20j], this implies that one can- 
not employ the solution of the moment problem ([5]) any 
longer. Here we will determine, on the basis of the 
Chebyshev-Markov and other model-free inequalities, the 
bounds for the asymptotic form of the fast projectile 
stopping power. Let us consider the contribution 

pk"<k'i jfc ra+ 

Si := / — / uj 2 nB (uj) C (k, uj) dui 



< j -^r«+ S(k,w)du, (7) 

on account of the fluctuation-dissipation theorem (FDT) 
17j . Then, by applying the upper bound obtained in 
Ref. [23| for the charge-charge static structure fac- 
tor of a quantal multi-component plasma under the 
assumption of perfect screening, limfc_o S(k)/k 2 < 
fkVp coih(hu> p f3 /2) / (8irne 2 ) , we can approximate the pre- 
vious integral as: Si < iyuj p coth(hujpf3/2)k"v, for k'{ < 
k[ ~ vf/v. 



This contribution should be compared with those stem- 
ming from 



S 2 := 
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— j lo 2 ub (uj) C(k,uj) duj 



uj C (k, uj) dui. 
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Clearly, we can find an upper bound for £2 analogous to 
expression ([6]). To determine a lower bound we might 
apply the Chebyshev-Markov inequalities (CMI) [l^]. In 
particular, if we take the measure da = uj 2 Cduj, then 



S 2 > 




(9) 



Since we have assumed that Vfc <E (k[, k%), kv > uj 2 (fc), 
then, for some £ > 1 such that t;k[ ~ Vf/v asv/vp ^> 00, 
we have 



S 2 > 




Hence, if we want to ensure that the first leading term 
of the stopping power asymptote is contained in £2 , it is 
sufficient to choose the lower cut-off as k'{ ~ v F /v 2 , to 
obtain Si < ttuj 2 coth(hujp/3/2)vF /v + O (y 2 F /v 2 ), which 
becomes negligible compared to S2 as v/vp — * 00. 

The last contribution to the stopping power ([2]) reads 



5, := 



k'J>k 2 



dk 



uj 2 ub (uj) C (k,uj) duj. (11) 



In particular, if k% ^> 2Mv/h = ^M/m, then 

X<S 3 < n B (.a-{k%))X, (12) 

X = dk/k uj 2 C (k,oj) dcu. By applying again the 
CMI, but now with the measure dS = Cdui, it is possi- 
ble to prove that the satisfaction of all three sum rules, 
{C2 V } 2 = , alone does not guarantee the convergence of 
53. To this aim we may introduce an additional condi- 
tion on the decay of the LF in the interval of interest 
(a_,aq-). Precisely, from the inequalities 



< / ui 2 dY l < a 2 



(13) 



we see that £3 converges if and only if dT, < (kp/k) 1 , 
7 > 4, which can be achieved by imposing on the distri- 
bution £ (uj) the following Holder condition: 



|£ («+)-£ («_)]< 





a + — a- 




UJ p 
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with < v < 1, /i > 3, for k > 

1 n - ' '■•>■"' Therefore, it is feasible to choose 



k%. Then, 



v 2 /vp to get that <S 3 < 
. A Holder-type condition 



S 3 < 2(v 2 p Mk 2 / {mk'{ 
an upper cut-off as k'2 ^ 
2uj 2 Mv f / (mv) + O (v 2 F /v 2 
like fT4")) can be fulfilled in a number of physical mod- 
els, namely: in an electron-ion hydrogen plasma, where 
the above-mentioned PE asymptote 18| is applicable if 
one assumes the spatial dispersion to be negligible for 
wavelengths much higher than the maximum impact pa- 
rameter, what is valid for the range of frequencies and 
wavenumbers considered for £3 . In case of a uniform elec- 



tron gas, the asymptotic expression derived in Ref. [24 1 
satifies a similar condition as well, although one needs 
to take into account the region of non-analyticity of the 
perturbative expansion (25| . 

With the aforementioned conditions, it follows that: 



Theorem 1 The stopping power 
satisfies asymptotically, as v/vp — 



"Vp 
V v 



< 



Vp 
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This theorem provides the bounds for the fast projec- 
tile asymptotic form leading term. These are based on 
inequalities which do not depend on the particular de- 
tails of the fluctuation spectrum at low and intermediate 
frequencies. 

Conclusions.- In this Letter we have studied the mod- 
ification of the BL expression for the plasma stopping 
power due to the presence of an ion component, strong 
coupling and the decay of the Langmuir mode. Wc 
have shown that, for a perfectly defined plasma collec- 
tive mode with negligible damping, the above-mentioned 
expression is affected by the electron-ion correlation. In 
addition, we have derived bounds for the fast projectile 
asymptotic, on the basis of well-established results of the 
linear response theory of Coulomb systems, namely, the 
zero-frequency sum rule, the f-sum rule, the fourth mo- 
ment sum rule, and the FDT, together with the com- 
pressibility sum rule. This general result constitutes a 
sum rule for the calculation or numerical estimate of 
the fast projectile stopping power for any model dielec- 
tric function satisfying the above-mentioned conditions, 
not only in plasma physics, but also in other multi- 
component uniform charged particle models of condensed 
matter physics. 
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